Generation of Noise Time Series with arbitrary Power Spectrum 
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Noise simulation is a very powerful tool in signal analysis helping to foresee the system performance 
in real experimental situations. Time series generation is however a hard challenge when a robust 
model of the noise sources is missing. We present here a simple computational technique which 
allows the generation of noise samples of fixed length, given a desired power spectrum. A few 
applications of the method are also discussed. 



INTRODUCTION 

A standard approach in noise time series simulation is 
based on Carson's theorem 1]. In fact, the superposition 
of randomly delayed pulses of a definite shape /(t) with 
arbitrary coefficients a*; gives rise to a pulse train 



n{t)^Y.akf{t-tk) 



(1) 



whose power spectrum Nioj) can be expressed in terms 
of the /(t) power spectrum Pf{i^) according to 



N{uj) = aPf{uj) = 



(2) 



where a is a proper normalizing constant and F{uj) 
is the Fourier transform of f(t). The delays t^'s are 
distributed according to Poisson statistics and f{t) is 
usually required to have a zero mean value so that n(t) 
averages to zero. It can be shown that this result is 
independent of the density function used for the random 
amplitudes a^, if the variance of such distribution is 
fixed by the choice of the normalizing constant a. 

The basic idea of the method consists then in searching 
a proper expression for f(t) so that the generated time 
series n{t) have the desired power spectrum N{uj). In 
other terms one is faced with the problem of inverting 
eq. ([2|) to obtain the expression of f{t) given Pf{Lo). 
Then, n{t) can be built iteratively according to eq. ([T]) 
with a proper choice of the random parameters Uk and tk ■ 

This approach can be quite difficult whenever the 
desired spectrum is different from a simple power law 
and analytical approaches to the determination of f(t) 
are usually not straight-forward. In fact, the desired 
noise spectrum often shows complex behaviours in 
the measured frequency range. An example of such a 
common situation is shown in fig. [1] (referring to the 
output of a few hundred microgram bolometric detector) 
where we have a complex forest of microphonic lines 
overimposed to a smoothly varying distribution. The 
anti-aliasing filter cut frequency is also apparent by the 
quick drop of the power spectrum in the 100-200 Hz 
region. 



The method proposed here consists in the determi- 
nation of an f(t) which exactly matches a given noise 
power spectrum. This is essentially based on the choice 
of a particular (out of an infinite number of possibilities) 
function f{t) which satisfies eq. ^ and can be applied 
both to the case in which N{uj) is known analytically or 
experimentally measured. 

Let us begin with the discrete case of a sampled se- 
quence /[fc]. If we consider a finite number L of samples 
and assume ergodicity, then the power spectrum can be 
approximated by an ensemble average over a sufficiently 
large number of finite sequences fL[k] of fixed length L 
according to 



Pf[k] ^< inhm'y 



(3) 



where J- stands for the Discrete Fourier Transform 
(DFT) operator. Since the sequence elements fL[k] are 
supposed to be real, the corresponding power spectrum 
is a real even function and only half of its values are in- 
dependent. In fact all the informations concerning the 
phases are lost in the quadrature and this is the reason 
why the inversion of eq. ([3]) does not admit a unique so- 
lution /i [fc] . 

The simplest solution to eq. ^ is obtained by deleting 
the ensemble average and inverting the resulting expres- 
sion. Phases can be then randomly added assuming a 
fiat distribution. The result represents the core of this 
work and can be summarized as follows: 



F[fc] = VN\k] ■ e'^" 



(4) 
(5) 



where iV[fc] is the power spectrum according to which 
we want to generate our time series (known a priori) and 
6k are random numbers uniformely distributed between 
and 27r. In order to guarantee the condition of reality 
of /[fc] upon inverse DFT the sYmmetry constraint 
F[k] = F*[—k] must be imposed. 

Noise time series can now be built following eq. ([l}, 
simply generating the delays tk according a Poisson dis- 
tribution: 



tk = tk-i - ln(l - R)/\ 



(6) 
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where R is a random number with a uniform distribution 
and A, which represents the overlap rate of the signals 
f[k], is the only free parameter of the method and is 
usually adjusted to improve the quality of the simulated 
time series. 

As stated above, the choice of the normalization term 
a partially fixes the arbitraryness in the choice of the 
amplitudes density function by setting its variance. This 
can be better understood by following some of the steps 
in the derivation of the Carson's theorem. Actually, by 
substituting eq. ^ into the the power spectrum defini- 
tion we obtain: 



N{io) = 



^akf{t - tk) 



(7) 



The last sum can then be separated into two different 
terms, one for i = j, the other for i ^ j. It can be shown 
that the latter is proportional to F(0) and therefore 
vanishes since we required /(t) to average to zero. We 
obtain therefore: 



iV(a;) 




(8) 
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FIG. 2: Superposition between a simulated (black) and mea- 
sured (red) jfet noise power spectrum. The simulated spec- 
trum was obtained averaging 50 finite time series. Error bars 
are the average standard deviations of the simulated spectra. 
The distribution of the x^/d.o.f. for a set of simulated spectra 
is shown in the inset. 
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FIG. 1: Superposition between a simulated (black) and mea- 
sured (red) bolometric noise power spectrum obtained with a 
bolometric detector of 0(mg) mass 5j. The simulated spec- 
trum was obtained averaging 50 finite time series. Error bars 
are the average standard deviations of the simulated spectra. 
The distribution of the x^/d.o.f. for a set of simulated spectra 
is shown in the inset. 



FIG. 3: Ratio between the simulated and measured bolomet- 
ric noise power spectra shown in Fig[T] 



N(L,) = l^F(L,)f^M <\af > (9) 

By comparing eq. ^ and eq. we can finally obtain 
the relation between the average rate A, the normaliza- 
tion constant a. and the amplitudes variance |ap : 



where T is the finite length of the time series. Any arbi- 
trary distribution with a variance given by eq. (jlOp can 
now be used to generate the random amplitudes Ofc and 
the method is fully defined. 

We can therefore summarize the most relevant steps 
for the method implementation as follows: 



The averaged sum on the right hand part can be 
thought as the average of the squared amplitude of 
each single pulse times their average number M of 
occurencies. As a result we can express the average 
power spectrum as: 



• Given the desired noise power spectrum N{uj) se- 
lect the basis function f(t) according to eqs. @ and 
([5]). Then eq. ([2]) fixes the constant a 
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number of averaged noise baselines 
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FIG. 4: Samples of simulated (top) and a real (bottom) noise 
time series of a jfet 
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FIG. 5: Ratio between a simulated and measured jfet noise 
power spectra shown in Fig[2] 



• Generate a set of increasing delays tk according to 
a Poisson distribution with average rate A. 

• Generate the random amplitudes oj. according to 
an arbitrary distribution with variance fixed by 
eq.dlOl). 

• Shift f{t) by tk imposing a periodicity constraint 
f{t) = f{t + T) and multiply the result by a^. 

• Sum iteratively over tk while ^ ifc < T 

The method can be used also when a theoretical 



FIG. 6: Average normalized standard deviation (dots) of the 
simulated power spectra as a function of the number of simu- 
lated time series used for the average process. The black line 
joins simulated points while the red one is an inverse square 
root law shown for reference. 



(analog) estimate of F{uj) is known. This is a common 
situation when the performance of an electronic circuit 
or of a detector are studied. In this cases the discrete 
power spectrum F[k] can be in fact obtained by means of 
a proper frequency warping 0] of the analog theoretical 
estimate of the power spectrum. 



APPLICATIONS 

We would like to discuss in this section a few exam- 
ples of application of the method to situations character- 
ized by a complex structure of the desired noise power 
spectrum. We have chosen to this end two completely 
different noise sources: a single particle bolometric de- 
tector and a jfet transistor Indeed, being 
characterized by very different masses, frequency ranges 
and noise sources, they provide a clear example of the 
potential and versatility of the method. 

The concept of a single particle bolometric detector 
is very simple. It is an ideal calorimeter for which all 
the energy of the incident radiation is converted to heat 
giving rise to a temperature rise of the detector's body. 
Such a temperature variation is then measured by means 
of a proper transducer which in our case will be rep- 
resented by a doped thermistor whose resistance varies 
with temperature. A bias circuit is therefore needed and 
the observed signals express the voltage change at the 
thermistor terminals. Working temperatures are as low 
as few tens of K in order to guarantee measurable tem- 
perature variations. In such working conditions differ- 
ent noise contributions can be recognized: the electronic 
noise coming from the acquisition chain, the Johnson 
noise in bias circuit elements, the mechanical vibrations 
or instabilities of the of the experimental setup which 
may produce variation in the acquired voltage samples 
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and finally the temperature instabilities due to thermo- 
dynamic fluctuations. The result is a complex noise 
power spectrum in which a forest of sharp lines is overim- 
posed to a generally smooth behaviour. Detector details 
(eg. mass) naturally affect the specific shape of the noise 
spectra (range, lines, ...) but usually do not alter the gen- 
eral trends described above [sj, |9| . Following the common 
prescriptions ^] such spectra were obtained by averaging 
a few hundreds of randomly selected baselines (defined 
as time windows in the absence of a signal). The same 
averaging procedure was followed to obtain a power spec- 
trum of the output noise of a simple jfet transistor (fig.[2]) . 
Following the recipe described in this paper and starting 
from the power spectra of a small mass bolometer and 
a jfet, a set of 50 simulated time series was generated. 
The respective power spectra were obtained according to 
eql3]by averaging their DFT's. The result is directly com- 
pared in fig. [1] and [2] with the original noise power spectra 
obtained experimentally. The ratio between the experi- 
mental and simulated power spectra is also reported in 
fig. [3] and in fig. [5l Experimental and simulated base- 
lines are also compared in fig. U in the case of the jfet 
transistor. 

In order to allow a more quantitative estimate of the 
level of agreement between the simulated and the origi- 
nal power spectra, a full set of simulated power spectra 
was generated in order to evaluate the distribution of the 
simulated values at each frequency. Each power spec- 
trum was obtained averaging 50 simulated time series. 
The obtained standard deviations are shown in fig. [2] and 



fig-Hi The distributions of the x^/d.o.f. are also shown to 
demonstrate the correct statistical behaviour. The fluc- 
tuations of the simulated power spectra values depend of 
course also on the number of time series used for the av- 
erage process (eqlH]). Such a dependence was explicitely 
studied by generating a different set of power spectra and 
evaluating the corresponding average standard deviation 
from the desired otput, as a function of the number of 
averaged time series. The result is shown in fig. [6] and 
found in excellent agreement with an inverse square root 
law. 
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